Majorana fermions emerging from magnetic nanoparticles on a superconductor 

without spin-orbit coupling 
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There exists a variety of proposals to transform a conventional s-wave superconductor into a 
topological superconductor, supporting Majorana fermion mid-gap states. A necessary ingredient 
of these proposals is strong spin-orbit coupling. Here we propose an alternative system consisting 
of a one-dimensional chain of magnetic nanoparticles on a superconducting substrate. No spin- 
orbit coupling in the superconductor is needed. We calculate the topological quantum number of a 
chain of finite length, including the competing effects of disorder in the orientation of the magnetic 
moments and in the hopping energies, to identify the transition into the topologically nontrivial 
state (with Majorana fermions at the end points of the chain). 

PACS numbers: 74.45.+C, 74.78.Na, 75.30.Hx, 75.75.Lf 
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I. INTRODUCTION 

The search for Majorana fermions in superconducting 
nanowircs is both rewarding and difficult These non- 
degenerate midgap states at the end points of the wire re- 
quire that the s-wave proximity effect coexists with bro- 
ken time-reversal and spin-rotation symmetries, which 
together can drive the superconductor into a topologi- 
cally nontrivial phased While the former symmetry can 
be readily broken by a magnetic field, it has been argued 
that Rashba spin-orbit coupling is too weak to effectively 
break the latter symmetry^ We are optimistic that some 
variation on the InAs nanowire proposal^ will be suc- 
cessfully realized, but alternative proposals^— continue 
to play an important role. 

In this paper we propose a route to Majorana fermions 
in s-wave superconductors that docs not at all require 
materials with spin-orbit coupling. We consider a one- 
dimensional chain of magnetic nanoparticles (magnetic 
dots-^) on a superconducting substrate. (See Fig. [1]) 
The nanoparticles create bound states in the supercon- 
ducting gap, having a nonzero magnetic momentJ^— 
The magnetic moment breaks time-reversal symmetry 
as well as spin-rotation symmetry, without the need for 
spin-orbit coupling in the superconductor. 

As we will show, the transition into the topologically 
nontrivial phase is governed by the competition of two 
types of disorder: On the one hand, disorder in the orien- 
tation of the magnetic moments on nearby nanoparticles 
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FIG. 1: Proposed setup to obtain Majorana fermions at the 
end points of a chain of magnetic nanoparticles (typically Fe 
or Ni) on an s-wave superconducting substrate (typically Pb). 
Arrows indicate the orientation of magnetic moments of the 
nanoparticles. 



is needed to open a gap in the excitation spectrum - 
which is a prerequisite for a topological phase. On the 
other hand, disorder in the hopping energies localizes the 
states on short segments of the chain — suppressing the 
superconducting order needed for the topological phase. 

In the next section we introduce the model of a chain 
of magnetic nanoparticles on a superconductor, and then 
in Sec. IIIII we use scattering theory^ to calculate the 
topological quantum number of a finite chain, including 
the competing effects of magnetic and hopping disorder. 
As shown in Sec. IIV1 we find a substantial region in pa- 
rameter space that is topologically nontrivial (and thus 
has Majorana fermions at the end points of the chain). 



II. MODEL OF MAGNETIC NANOPARTICLES 
ON A SUPERCONDUCTOR 



A. Hamiltonian 

For each magnetic nanoparticle on top of the super- 
conductor we consider a single electronic orbital near the 
Fermi level fi, with spin a coupled to the local magnetic 
moment through an effective magnetic or exchange field 
B n . We set h and the Zccman energy gfiB both equal 
to unity, so B n is measured in units of energy. We as- 
sume that the magnetic moment has the same magnitude 
Bo = \B n \ on each nanoparticle, varying only in the ori- 
entation. States on neighboring nanoparticles n, n+ 1 arc 
coupled by a hopping energy t„ , where the n-dependence 
accounts for variations in the nanoparticle spacing. 

The superconducting substrate induces a spin-singlet 
pairing energy Ao on the nanoparticles, taken as n- 
independent. The electrostatic potential induced by the 
superconductor is also taken as n-independent, so that it 
simply gives an offset to fi. The charging energy e 2 /C is 
unimportant because the nanoparticle is strongly coupled 
to the superconductor. 
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The Hamiltonian of this model is given by 

H = (Wna/n+l.a + H - c ) ~ ^Y fnafna 

+ ^ (B n ■ (T) a pfl a f n i3 

+E( A °44+ Rc -)' (2 - x) 



where the abbreviation H.c. stands for Hermitian conju- 
gate. The operator f na is the fermion operator for a spin- 
a electron on the n-th nanoparticle and cr = (a x ,a y ,a z ) 
denotes the vector of Pauli matrices. The magnetic mo- 
ments are taken as frozen, without any dynamics of their 
own. 

Upon transformation to the Bogoliubov basis ^ n = 
(/nt) fni, flp ~f}v\)i trie Hamiltonian takes the form 



H = H [{HtnVn+i + H.c.) + *t h n V n 



(2.2) 



The matrices h n and i n are defined by 
h n = 



-fxa + B n ■ cr A a 

A <t A^o + B n er 
LivoTz + {B n ■ <t)t + A a T x . 











= t n a T z 



(2.3) 
(2.4) 



The Pauli matrices <jj and Tj (i = x,y,z), with identity 
matrices <jq and To, act on the spin space and particle- 
hole space, respectively. 

The Hamiltonian (|2.2[) has neither time-reversal nor 
spin-rotation symmetry, but it does satisfy the particle- 
hole symmetry requirement that 



H I y -H fOr H-» (JyTy^. 



(2.5) 



This places the system in symmetry class D, which in one 
dimension has a topologically nontrivial phased 



with coefficients defined in App. [A] The effective pair 
potential A n is of order A t n /B , dependent on the rel- 
ative angle between B n and B n +\. For parallel magnetic 
moments A n vanishes. 

The single-band Hamiltonian (|2.6[) has the same form 
as Kitaev's spinless p-wave superconducting chain^ The 
difference is that here the p-wave pairing is obtained from 
s-wave pairing due to the coupling of the electron spin 
to local magnetic moments. This has the effect of cou- 
pling spin to orbital degrees of freedom, but in contrast 
to existing proposals^— neither spin-orbit coupling in 
the superconductor is needed nor a Rashba effect in the 
nanowire. 

We have given the single-band limit (|2 to make con- 
tact with Kitaev's model. In what follows we will use the 
full Hamiltonian (|2.2[) . valid to all orders in 1/Bq. 



C. Disorder 

We distinguish the localizing effect of disorder in the 
hopping energy £„, which localizes the electrons without 
opening an excitation gap, from the gap opening effect of 
disorder in the orientation b n = B n /Bo of the magnetic 
moments. 

Disorder in the hopping energies, due to variations in 
the interparticle spacing, is modeled by drawing the i n 's 
(n = 1, 2, . . . N — 1) independently from a uniform distri- 
bution in the interval (to — Stto,to + 5tto). 

For the magnetic moments we take a dipolar ferromag- 
netic correlation of the unit vectors b n on neighboring 
nanoparticles, according to the distribution 

P{bi,b 2 ,...b N ) = (4tt) II , , 

n=i Sbsmh6 b 

(2.7) 

where the parameter 5b > quantifies the strength of the 
correlation (strongly correlated for small <5&). 



III. SCATTERING MATRIX 



B. Single- band limit 

In the large magnetization regime Bq ~ 3> Ao,t n 
the electron spin on the n-th nanoparticle is nearly po- 
larized along B n . As outlined in App. [A] the Hamilto- 
nian can then be projected onto the lowest spin band, 
with electron operator ip n . Only virtual transitions to 
the higher spin band contribute. To first order in I/Bq 
the effective Hamiltonian has the form 



H eS = y 



(t n 1p^1p n+ l + t^VnVWa + H.C.) - jln^n^n 



(A„vM+i+H. 



(2.6) 



To identify the topologically nontrivial phase of a finite 
disordered chain of N nanoparticles it is more efficient to 
work with the scattering matrix than the Hamiltonian J£ 
The scattering matrix <S relates incoming and outgoing 
wave amplitudes at the Fermi level. The waves can come 
in from the left end or from the right end of the chain, 
in two spin directions and as electron or hole, so S is 
an 8 x 8 unitary matrix. Its 4 x 4 sub-blocks are the 
reflection and transmission matrices, 



S = 



71 V 
T TV 



(3.1) 



Particle-hole symmetry (|2.5p requires that 

X = a y T y X*a y T y , X e {71, 71' , T, T'}. (3.2) 
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Following Rcf. |20j we calculate the scattering matrix 
by writing the tight-binding equations at the Fermi level 
in the form 



t n <f>n\ _ /i„-l$ n -l 



M, 



o t n 



(3.3a) 
(3.3b) 



where $ n is a four-component vector of wave amplitudes 
on site n. (Sites n = and n = N + 1 represent electron 
reservoirs.) Waves at the two ends of the chain are related 
by the transfer matri»2i 

M = M N M N -i---M 2 Mi. (3.4) 
Current conservation is expressed by the identity 

'0 -i" 



i 



(3.5) 



We transform to a new basis with right-moving and 
left-moving waves separated in the upper and lower four 
components, by means of the unitary transformation 



1 A 1 

i —i 



(3.6) 



The current conservation relation (13.51) transforms to 



MtX z M n = X z , E 2 



1 

-1 



(3.7) 



In this basis the transmission and reflection matrices fol- 
low from 



(3.8) 



Unitarity together with particle-hole symmetry 
ensure that the determinants of 71 and 71' are identical 
real numbers. The Z2 topological quantum number Q = 
±1 is given b y 18 ' 22 



Q = sign Det ft. 



(3.9) 



As shown in Ref. [01 Majorana bound states exist at the 
end points of the chain if and only if Q = — 1, so this 
identifies the topologically nontrivial phase. 



IV. PHASE DIAGRAM 




FIG. 2: Determinant of the reflection matrix as a function 
of n, for fixed Bo = 2to and A = 0.9 to, ensemble averaged 
over 400 chains of nanoparticles of length N. This is data for 
St — 00, St — 0, corresponding to random and uncorrelated 
orientations of the magnetic moments of the nanoparticles, 
without randomness in the hopping energies. 



realization and a single system size (N = 6000). The N- 
dependence of the ensemble average (Det r) is plotted in 
Fig. [2] to show how the transition becomes sharper with 
increasing system size. 

The competing effects of the two types of disorder are 
evident in the phase diagrams. We see that the topo- 
logically nontrivial region is largest for <5 t equal to zero 
and 5b small but not too small. Disorder in the hopping 
energies or in the magnetic moments causes localization, 
which reduces the nontrivial region (mainly on the small 
Ao side). Some disorder in the orientation of the magne- 
tization is needed to open a gap in the spectrum, so while 
the optimal value of S t equals zero, the optimal value of 
St is small but nonzero. 

To obtain more insight in the Ao and /i dependence of 
the phase diagram, we have calculated the phase bound- 
ary by the method of self-consistent Born approxima- 
tion (scba). As shown in Ref. [23| for a different system, 
the SCBA can locate the gap inversion associated with 
the topologically nontrivial phase. The calculation is de- 
scribed in App. [B] and the result is shown in Fig. [3] (red 
dashed curve in top panel) , for a system without disorder 
in the hopping energies (5t = 0) and with random and un- 
correlated directions of the magnetic moments (5b = 00). 

The SCBA describes reasonably well, without any ad- 
justable parameters, the location of the phase boundary 
at large A . The phase boundary at small A is not re- 
covered by the SCBA. Our explanation is that the former 
phase boundary is due to gap inversion, while the latter 
phase boundary is due to localization (which is beyond 
the reach of the SCBA). 



Results of our search for a topologically nontrivial 
phase are shown in Figs. [5]-[3] We have kept B n = 2t 
fixed and varied the material parameters Ao , /i as well 
as the disorder parameters 5b, 5t- The color scale Figs. 
|3] and 2] of Det r show the topologically nontrivial region 
emerging as a black island (Det r = — 1) in a yellow back- 
ground (Det r = 1). These are plots for a single disorder 



V. CONCLUSION 

In conclusion, we have proposed to create Majorana 
fermions by depositing magnetic nanoparticles on a su- 
perconducting substrate. The superconductor has con- 
ventional s-wavc pairing and need not have any spin- 
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FIG. 3: Color scale plot of the determinant of the reflection 
matrix as a function of [i and Ao, for fixed Bo = 2to, calcu- 
lated for a single disorder realization in a chain of N = 6000 
nanoparticles. The topologically nontrivial phase has Detr = 
— 1 (black region). The four panels are for different values of 
the disorder St in the hopping energies and Sb in the orien- 
tation of the magnetic moments [with Sb = oo corresponding 
to random and uncorrelated orientations, see Eq. (|2.7|l ]. The 
dashed red curve in the top panel gives the phase boundary 
following from the self-consistent Born approximation. 



orbit coupling. The candidate material for the magnetic 
nanoparticles is iron or nickle. For the superconductor 
one could choose lead, with a coherence length £ ~ 80 ran. 

The magnetic moments of nearby nanoparticles (spac- 
ing small compared to £) are correlated by dipolar 
interactions^ As we have found, a strong correlation in 
the orientation of the magnetic moments (Sb <C 1) does 
not preclude the appearance of a topologically nontriv- 
ial phase — only a nearly complete alignment suppresses 
it. The most stringent requirement seems to be the need 
to have little disorder in the hopping energies, since in 
this one-dimensional system the localization length is of 
the order of the mean free path and hence severely lim- 
ited by disorder. The spacing of the nanoparticles should 
therefore be highly uniform. 

While a long chain of nanoparticles (much longer than 
£) is needed to create a pair of weakly coupled Majorana 
fermions at the end points, a signature of a topological 
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FIG. 4: Same as Fig. [3] but now as a function of Sb and St for 
fixed Ao, At. 




FIG. 5: Excitation spectrum (black solid curves) of two mag- 
netic particles with an angle 9 — 70° between their mag- 
netic moments, calculated from the Hamiltonian (|2.1[l at 
fixed fi = Bo = 2Ao as a function of the hopping energy 
to- The level crossings at the Fermi energy do not split be- 
cause the ground states at the two sides of the crossing dif- 
fer in fermion parity V SE {0, 1} [red dashed curve, calcu- 
lated from the Pfaffian of the antisymmetrized Hamiltonian, 

p = i-i s igiiPfK7vff)]. 



phase transition can be identified even if one has only 
a few nanoparticles. This signature is a protected level 
crossing at the Fermi level, which signals a change in the 
fermion parity of the ground state.— We show this feature 
in Fig. for the simplest case of just two nanoparti- 
cles with a non-collincar magnetization. Nanoparticle- 
induced bound states within the superconducting gap 
have been resolved by scanning tunneling microscopy, 
for a Mn dimer adsorbed on a Pb superconducting thin 
filmj2£ We predict level crossings at the Fermi energy 
(E = 0) upon variation of the spacing of the two mag- 
netic atoms, as a signature of the switch from an even to 
an odd number of electrons in the ground state. 
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Appendix A: Relationship to Kitaev's model 

We establish the relationship between the tight- 
binding model (|2.1j) of electrons coupled to an s-wave 
superconducting order parameter and to a randomly ori- 
ented magnetization, on the one hand, and Kitaev's 
model of a spinless p-wave chain^ on the other hand. 

The magnetic moment B n = B b n on the n-th 
nanoparticle has fixed magnitude B and orientation 
b n = (sin#„ cos (f) n , sin(9„ sin </>„, cos 9 n ). We align the 
spin basis on each site with b„, 



fir] 



Un 



_ ( cos(#„/2) -sin(0 n /2)e-^™ 
sm{e n /2)e 1 ^ cos(0 n /2) 



The Hamiltonian (12.11) transforms to 



(Al) 
(A2) 



+ ^2 B a z.,a{3friafn{3 ~ fnafna 
n,a,/3 n.a 

+E( A °/^4+ H - c -)- ( A3 ) 

n 

The unitary matrix f2„ = U^U n +i has elements 



a n = cos — cos 



a n -ft* 

fin Oi n 



— rns V- sm — sin p vr 



(A4a) 



2 2 



sm — sm — - — e 

2 2 



B„ = - sm — cos — — e 9n + cos — sm — 
2 2 2 2 



(A4b) 
+i _ 
(A4c) 



In the regime B — \n\ 3> t n , A, only one of the two 
spin bands lies near the Fermi level. (For definiteness, we 
take /x < 0, so that the spin-down band is the band near 
the Fermi level.) The effective low-energy Hamiltonian 
-ffeff has only virtual spin-flip transitions, without any 
matrix elements between the spin bands. We obtain H g 
from H by a canonical transformation H ^ e~ lS He lS 
followed by a projection onto the spin-down bandi2S~— 
The Hermitian operator S is expanded as S = Si + S2 + 
• • • , with S n of order A^Bq T n~ s and p + q = r + s = n. 
The term S n is chosen such that it eliminates the n-th 
order matrix elements between the spin bands. 



H ^ H + [H,iSi], (A5) 



2^/ni/»t-H.c. 



The resulting effective single-band Hamiltonian is 

tn\ 2 + 



E ( BtH 

n 

E 



2B 



, * ,t / 1 tn+lfin+ltnftn , + , 



(A6) 



n+2 



1 

2B~ 



1 

2fi 



Aotufinlpnlp 



n+1 



H.c. 



, (A7) 



where we have abbreviated rp n = f n ^. This is the p- 
wave chain model of Kitaev^ with the addition of a next- 
nearest-neighbor hopping term. The pair potential cx 
Aoi n /3 n vanishes if the magnetic moments are aligned, 
because then /?„ = 0. 



Appendix B: Phase boundary in self-consistent Born 
approximation 

We calculate the phase boundary of the topologically 
nontrivial phase using the self-consistent Born approxi- 
mation (scba), starting from the Hamiltonian (|A3|) in 
the locally rotated spin basis. For simplicity we take 
site-independent hopping energies to and random, un- 
correlated magnetic moments. 

On each site n we define four real random variables 

(n) , v? 2 H , ¥>3 (™) , ¥>4 («■) 

= Rea„ - |, Im.a n ,Refi n ,Im/3 n , (Bl) 

of zero average and variance 

(<Pj(n)<pj>(ri)) = Sjj'Snn'-^-, 



(TO 1 ,m 2 ,TO3,m 4 ) = (|,8,4,4). 



(B2) 
(B3) 



We neglect higher order cumulants, approximating the 
distribution of the ipj (n)'s by a Gaussian. (These random 
variables can be thought of as bosonic fields of mass rrij.) 

We work in the Bogoliubov basis vj/ n = 
(fnt,fni,.fl v -f nt ) and transform variables from 
site indices n G {1,2, ...iV} to dimcnsionless momenta 
k G (— 7T,7r). Substitution of Eq. (|A3I) gives 

4 

H = ^*+if (fc)* fe + ^^^(fc-g)*K(fc,g)* g , 

k j=l k,q 

(B4) 

Ho(k) = (§io cosk - p)aoT z + B a a z T + A a T x . (B5) 
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Wc have abbreviated 



N 



k 11 



(B6) 



The interaction vertices are given by 

v J (k,q) = ±t Q (a J e lq + a]e- lk ), (B7) 

«i = &ot z , a 2 = icr z T z , a 3 = -i<J y T z , a 4 = o- x t . 

(B8) 

They satisfy Vj (k, q) — Vj (q, k) and the particle-hole sym- 
metry condition 



j(k,q)a y 



The self-energy at the Fermi level is given in SCBA 
by the integral equation 



(k, q )[H (q) + nq)r 1 v J (q,k). 



(BIO) 
(Bll) 



The effective Hamiltonian 

H SCBA {k) = H (k) +23(fc) 
satisfies the particle-hole symmetry relation 

<TyT y H SCBA (k)(T y Ty = -ff s * CBA (-fc). (B12) 



In contrast to the SCBA calculations in Refs. I23ll29l . 
where the interaction vertices are momentum indepen- 
dent, the randomness of the magnetic moments induces 
momentum-dependent vertices in Eq. (|B4|) . As a conse- 
quence the self-energy becomes momentum depen- 
dent. 



We calculate the topological quantum number Q from 
-ffscBA by means of Kitaev's Pfafhan formula,— 

Q = sign {Pf [a y TyH SCBA (0)] Pf [a v TyH SCBA (ir)]}. 

(B13) 

The particle-hole symmetry relation (|B12[) ensures that 
the Pfaffian is calculated of an antisymmetric matrix. 
In the large- N limit this Hamiltonian expression for the 
topological quantum number is equivalent to the scatter- 
ing matrix expression p.9p ^ 

The 4x4 self-energy matrix can be decomposed in 
terms of 16 combinations of Pauli matrices (TqT^, 



(B9) £(fc) = Si(fc)CT T z + E2(fc)CT z To+E 3 (fc)(7oT, + - • • . (B14) 



At k = 0,7r only the terms shown are nonzero. These 
three terms appear in H SCBA (k) as a renormalization of 
the energy scales 

p,(k) =n-Hi(k), B(k) = B + Z 2 (k), 

A(k) = A + S 3 (fc). (B15) 

Substitution into Eq. (|Bl~3l) gives the topological quan- 
tum number 

Q = sign 



{ [B(0) 2 - e(0) 2 - A(0) 
B(tt) 2 - e» 2 - A(^) 2 ] } , 



e(k) = |fo cos k — p.(k). 



(B16) 
(B17) 



The phase boundary where Q changes sign is indicated 
in Fig. |3] (top panel). This calculation contains no ad- 
justable parameters. 
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